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Abstract 

In previous work [L. Blanchet and A. Le Liec, Phys. Rev. D 80, 023524 (2009)], motivated by 
the phenomenology of dark matter at galactic scales, a model of dipolar dark matter (DDM) was 
introduced. At linear order in cosmological perturbations, the dynamics of the DDM was shown 
to be identical to that of standard cold dark matter (CDM). In this paper, the DDM model is 
investigated at second order in cosmological perturbation theory. We find that the internal energy 
of the DDM fluid modifies the curvature perturbation generated by CDM with a term quadratic 
in the dipole field. Lhis correction induces a new type of non-Gaussianity in the bispectrum of 
the curvature perturbation with respect to standard CDM. Leaving unspecified the primordial 
amplitude of the dipole field, which could in principle be determined by a more fundamental 
description of DDM, we find that, in contrast with usual models of primordial non-Gaussianities, 
the non-Gaussianity induced by DDM increases with time after the radiation-matter equality on 
super-Hubble scales. This distinctive feature of the DDM model, as compared with standard 
CDM, could thus provide a specific signature in the CMB and large-scale structure probes of 
non-Gaussianity. 
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I. INTRODUCTION 



The concordance cosmological model ACDM, based on a cosmological constant A and 
cold dark matter (CDM), is the state-of-the-art of today's cosmology [lj. One of the most im- 
portant successes of ACDM is its ability to reproduce the observed acoustic power spectrum 
of temperature fluctuations in the cosmic microwave background (CMB) @, Q. However, 
despite many successes at cosmological scales, the ACDM model encounters some problems 
at galactic scales (down to a few fcpc), notably related to rotation curves of galaxies and the 
Tully-Fisher relation (see Refs. (Ulsl for reviews). In particular, it has been shown that the 
challenging observations for CDM at galactic scales can be summarized in a single empirical 
law, still unexplained today, and dubbed the modified Newtonian dynamics (MOND) 

The MOND proposal fostered many attempts at finding a relativistic modification of 
gravity, i.e. an extension of general relativity (see e.g. Refs. |9l-ll2l|). Alternatively, based on 
an analogy between MOND and the physics of dielectric media, the concept of dipolar dark 
matter (DDM) was introduced, at the Newtonian level, as a reinterpretation of MOND in 
terms of polarizable dark matter particles without a modification of gravity 13| . A relativis- 
tic model of DDM in standard general relativity was then proposed [l4 , 15|| . In this model the 
dark matter particles are endowed with a space-like vector field called the "dipole" moment. 
It was shown that the DDM has the potential to reproduce the phenomenology of MOND 
at galactic scales, and is strictly equivalent to ACDM at the level of first-order cosmological 
perturbations around a Friedman-Lemaitre-Robertson- Walker (FLRW) background. Thus, 
the DDM behaves like ordinary CDM at early cosmological times. Furthermore the model 
naturally involves a cosmological constant. 

In the present paper we study the DDM model 14, 15[ up to second order in cosmological 
perturbation and show that, at this order, the non-linear dynamics of DDM starts departing 
from that of ordinary CDM. Consequently, we find that DDM predicts an additional contri- 
bution to the curvature perturbation, specifically given by the internal energy of the DDM 
fluid. We analyse this extra contribution in the evolution of the curvature perturbation at 
second order in perturbation theory within the Langlois-Vernizzi formalism (l6l . 171 . 



Non-Gaussianity of primordial perturbations in general (multi-field) models of inflation 



has received a lot of attention in the last few years (see e.g. Refs. [18|, [19| for reviews). The 
best current constraints on non-Gaussianity in the CMB are from the WMAP7 data with the 
interval — 10 < /nl^ 1 < +74 (95% CL) on the usual parameter J^l 3,1 that characterizes the 
amplitude of the simplest type of non-Gaussianity, namely the local shape. This constraint 
should be improved by the Planck satellite, and may exclude single-field inflation models if 
a nonzero value is favored by the data. 

Here we point out that alternative models for dark matter, such as DDM, can provide 
another possible source of non-Gaussianity; in the case of DDM because the dipole field 
yields a second-order contribution to the curvature perturbation. We compute the bispec- 
trum of the curvature perturbation, limiting ourselves to super-Hubble scales, and find a 
specific contribution to non-Gaussianity due to the DDM, assuming that the initial statisti- 
cal distribution of the dark matter dipole moments is Gaussian at early times. We use two 
methods, one based on the one-loop bispectrum and the vanishing of the dipole moment 
in the background (that assumption being consistent with our previous perturbative de- 
scription of DDM); the other method based on the tree-level bispectrum and assuming that 
there is a small homogeneous background value for the dipole moment. The two methods 
are found to agree on the result when interpreted in the sense of cosmic variance. 
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The amplitude of the DDM non-Gaussianity signal depends on the a priori unknown value 
of the dipole moment at early times. However, we find that in contrast with usual models of 
primordial non-Gaussianities, where the curvature perturbation is frozen on super-Hubble 
scales during the standard cosmological era, the amplitude of the non-Gaussianity induced 
by DDM increases with time after the radiation-matter equality, because the dark matter 
dipole field evolves after this epoch. Although a full study of the evolution of perturbations 
on sub-Hubble scales would be necessary for detailed predictions, we expect this new type 
of non-Gaussianity to lead to a very specific signature, which could be identified by using 
different cosmological probes such as the CMB and large-scale structure surveys. 

The plan of this paper is the following. In the next Sec. [Til we present the model of dipolar 
dark matter. In the subsequent Sec. IHH we give a perturbative description of the dipolar 
dark matter up to second order in cosmological perturbations. In Sec. IIV| we compute the 
evolution of the curvature perturbation up to second order using a standard perturbation 
around a FLRW background. Finally, we compute and analyze the non-Gaussianity induced 
by the DDM contribution to the curvature perturbation in Sec. |V] Throughout this paper 
we used "geometrized units" where G = c = 1. 



II. MODEL OF DIPOLAR DARK MATTER 

Thephenomenological model of dipolar dark matter (and dark energy) developed in 



Refs. [14], [15[ is based on a matter action in standard general relativity, 

Sbdm = J d 4 x L DDM [ 9lMU , J", e] , (2.1) 

which is to be added to the Einstein-Hilbert action of general relativity (without a cosmolog- 
ical constant), and to the actions of all the other matter fields and fluids (such as baryons, 
photons, neutrinos, etc.), all described in the standard way. The DDM fluid in a space-time 
with metric g^ u is described by: 

1. A conserved mass current J M = au^, such that V^J^ 1 = 0, where u p is the four- velocity 

1 /2 

normalized according to u^u^ = —1, and a = (— J M J M ) is the rest mass density; 

2. A dipole moment vector field £ M , which intervenes in the dynamics only through its 
projection orthogonal to the four- velocity, £^ = L^^, where _L MI/ = g^ u + u^u v . From 
this vector we construct the polarization H"^ = cr£^, i.e. the number density of dipoles, 

and its norm n ± = a£ ± , with £ ± = (±^^) 1/2 = (g^&Q) 1 ' 2 . 



The Lagrangian describing the dipolar dark matter fluid reads [141 . |15 



J DDM 



-a + J,e-W(U ± ), (2.2) 



where the dot stands for the covariant derivative with respect to proper time, namely £ M = 
u v V v ^. Notice that we have = J^ ± + V p (J p u v £ u ), so that the Lagrangian fl2T2|) 

admits an equivalent form depending only on the orthogonal projection £jj this shows that 
the only dynamical degrees of freedom of the dipole moment £ M are those of the projection 
£j\ which is a space-like vector. Notice also that the first term in Eq. ( 12. 2 p is the Lagrangian 
of a pressureless perfect fluid, i.e. that of ordinary CDM. 
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The potential "W is a function of the norm n± = (T^ITTF) 1/2 of the polarization field. 
Its expansion in powers of U± is determined, in the weak-field limit H± <C ao only, by the 
requirement of recovering the phenomenology of MOND in the non-relativistic regime. Up 
to third order inclusively, it reads 

A 1 flTT 2 

5T(n ± ) = — + 2vr ni + — ni + o(u\) , (2.3) 

07T ~~ oao 

where A is the cosmological constant, and ao is the constant MOND acceleration scale, 1 
supposed to be universal and measured to the value ao — 1.2 X 10~ 10 m • s~ 2 0,01. In the 
strong- field regime n_i_ 3> a (but still non-relativistic), the potential W can be adjusted so 
as to recover the ordinary Newtonian limit (lij . 

Varying the Lagrangian (12. 2p with respect to the dipole moment one obtains some 
non-geodesic equations of motion for the dipolar fluid: 

u» = _|£ W' , (2.4) 
with the notations W = d^/dlli. and £^ = £j_/£±. Varying with respect to the conserved 



current is achieved using a convective variational approach (see the Appendix of Ref. [15]) . 
One obtains, after defining f2 M = ^ + u M (1 + 2£±W), 

ty, = -V^OT - n ± 5T') - R WVT uT±u T , (2.5) 

O" 

where the second term in the right-hand-side (RHS), which involves the Riemann tensor, is 
analogous to the standard coupling to curvature for the motion of particles with spins in an 
arbitrary curved background. Finally, the stress-energy tensor of the dipolar fluid can be 
derived by varying the Lagrangian ( 12. 2ft with respect to the metric, keeping = ^/—gJ^ 
as independent of the metric (since J£ is conserved in the ordinary sense, i.e. d^J^ = 0). 
One obtains 

T nv = Q ( PJ u) _ v p ( Ii p L u^ - u p Ii { l u u) ^j -(W- H ± W') g pv . (2.6) 

This stress-energy tensor is conserved, VjyT Miy = 0, as a consequence of the equations of 
motion (12 ,4p and evolution ( I2.5p . 

Up to an hypothesis of "weak clusterisation" of dipolar dark matter during the cosmolog- 
ical evolution [14|, it was shown that the model reproduces the phenomenology of MOND 
in the non-relativistic regime, in the sense that the Bekenstein-Milgrom [20] modification of 
the Poisson equation is recovered. The Euclidean norm of the ordinary gravitational field 
is then given by the derivative of the potential, g = W(U.±), and the MOND interpolating 
function is related to the potential by /x = 1 — 4n^(g)/g, where $(g) is the inverse function 
of W{U ± ), i.e. is such that U ± = $(#). 



1 Appearing both in the expansion (12.31) . the constants A and ao should be related numerically in this 
model, i.e. A ~ a§ (in good agreement with observations). 
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III. PERTURB ATIVE DESCRIPTION UP TO SECOND ORDER 



A. Equivalence between DDM and ACDM at first order 

In the context of cosmological perturbations around a FLRW background, the dipole 
moment, which is a space- like vector field, has to be perturbative in order not to break the 
isotropy of the background metric. 2 We thus write £^ = 0(1), which from the equation 
of motion (j2.4j) . together with Eq. ( 12. 3p . gives ix^ = 0(1). In the following, we will write 
apart the term corresponding to the minimum of the potential (12. 3p and playing the role 
of the cosmological constant (dark energy). Accordingly, we split the stress-energy tensor 
(I2.6P into T^ 11 = — ^ g^ u + T^ M , and focus on the dipolar dark matter contribution Tj^ M . 
It was shown in Ref. [14| that it is possible, up to first order in cosmological perturbations, 
to recast the stress-energy tensor of this fluid in the form of that of a pressureless perfect 
fluid, by means of a perturbative redefinition of the four-velocity. Indeed, we define as a 
convenient variable (Jz^ denotes the Lie derivative along m m ) 

K » = # l & = £l-gy v u», (3.1) 

which is 0(1) and satisfies k^u^ = —^±Ufi = 0(2). Note that Eq. (13. ip is intended to be (by 
definition) valid at any perturbation order. Then, we make use of the following redefinitions 
of the energy density and four-velocity at first order: 

p = a- V U U1 + 0(2) , (3.2a) 
If = u» + « M + 0(2) . (3.2b) 

Although the new redefined energy density p includes a dipolar contribution, it represents 
more importantly the conserved energy density associated with the new redefined velocity 
field U^, i.e. V^J7" M = 0(2) where = pU^ is the new current. One can then check 
that the stress-energy tensor of dipolar dark matter becomes identical to that of a cold dark 
matter fluid: 

T&u = P WW + °( 2 ) = T cdm + 0{2) . (3.3) 

Therefore, at the level of first-order perturbations, the dipolar fluid (which also involves the 
cosmological constant) is equivalent to the concordance model ACDM. 

This equivalence between DDM and ACDM can also be proved directly at the level of 
the Lagrangian ( 12. 2p . Indeed, at first perturbative order, the redefined conserved current 
Jp = pU^ is related to the original one, J M = au^, by 

= J» + atf - u^VuIil + 0(2) . (3.4) 

______ ^2 

Computing the norm of the current (13.4)) . say r = (—J^J 11 ) , and using the normalization 
condition = —1, we find 

r = a-J,e ± - + 0(2) . (3.5) 



2 In one approach to computing the non-Gaussianities, in Sec. IV CI below, we shall relax this hypothesis 
and assume a small homogeneous background value for the dipole moment. 
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We use this result into Eq. (12. 2p and discard the covariant divergence V^LT^ since it does 
not affect the dynamics. 3 Finally, neglecting also the terms quadratic and cubic in the 
polarization II j_ = er£j_ in the potential W, we can recast, at first order in the dipole 
moment, the Lagrangian (j2.2p - (l2.3p of our DDM model as 

£ddm = -r - — + 0(2) = L ACDM + 0(2) , (3.6) 

87T 

and we see that it agrees with the Lagrangian describing the dark sector of the cosmological 
model ACDM — namely a mass term r without interaction and the cosmological constant 
(neglecting other species than dark matter and dark energy). 



B. Energy density, pressure, and anisotropic stresses at second order 

We now extend the previous analysis to second order in cosmological perturbations, 
where, as we shall see, deviations from CDM start occurring. The explicit form of the 
quadratic term in the expansion ( 12. 3 \ of the potential yields H±_W = inYl^ + 0(3), such 
that W - U ± W = ± - 2tt U 2 ± + 0(3). From the requirements that U^U» = -1 + 0(3) and 
that the current = pU^ be conserved at this order, V ^J^ = 0(3), we find the following 
second-order accurate expressions [generalizing Eqs. (13. 2p ] for the new variables p and 
as functions of the initial variables a, u^, and 

p = a -\/ u U u ± -AnU 2 ± -^K 2 + 0(3), (3.7a) 

I/m = (i + ^ni + !j\ u»+ (l + k? + 0(3) . (3.7b) 

Here where k v is defined by Eq. ( 13. ip . valid to all orders in perturbation theory. 

We can then perform a standard 3 + 1 decomposition of the stress-energy tensor with respect 
to the redefined four- velocity U^, and obtain 

t ddm = (e + p ) U^U V + Pg» v + W v + 0(3) , (3.8) 

with the following second-order accurate expressions for the energy density e, pressure P, 
and anisotropic stress tensor jy iv : 

e = a- V U U U ± - 2vr - ok 2 , (3.9a) 

p = f ni - , (3.9b) 

TT = -4tt U^U^ - gk^k u) . (3.9c) 

The brackets ( ) stand for the symmetric and trace-free (STF) part of a given spatial tensor 
(orthogonal to C/^), e.g. k^k^ = k^k u — |k 2 _L^; from now on = g^ v + U fl U u is defined 



3 Notice that the second term in the RHS of Eq. (j3.5j) can easily be transformed, up to a total divergence, 
into crUfj,^. From the equations of motion w M — 0(1), this term would seem to be of second perturbative 
order. However, one is not allowed to use the equations of motion back into the Lagrangian, so the second 
term in (13.51) is to be viewed as first order when considered in the DDM Lagrangian (12.21) . 
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with respect to the four- velocity ( 13.7bj) . The anisotropic stress tensor 'E fJ,u is orthogonal 



to the four-velo city U » and traceless, i.e. U U W U = 0(3) and = 0(3). Notice that 
the redefinition ( I3.7bj) of the four-velocity has allowed us not only to define the associated 
conserved mass density (I3.7a|) . but also to cancel any possible second-order "heat flow" Q M : 
there is no term of the type 2Q^U V \ with Q^W = 0, in Eq. (ETS|) . 

Now, we can introduce the "specific internal energy" W of the fluid via the definition 
e = p(l + W). This dimensionless quantity will play an important role in the following; it 
explicitly reads 

W = —U 2 ± -^- + 0(3). (3.10) 
o I 

Hereafter we substitute the four- velocity = w M + 0(1) in place of in the definition 
(13. ip of since the difference will only appear at third order in perturbation theory. The 
projected conservation of the stress-energy tensor, U fJj 'V u T^ )M = 0, yields 

e + Q(e + P + f3) = 0(3), (3.11) 

where G = V M £/ M is the expansion of the fluid, and = E^o^/O the so-called dissipative 
pressure, with a^ u = V ' (nU v \ the STF shear tensor. The dot now stands for the covariant 
proper time derivative along the new fluid worldlines with tangent four- velocity If 1 , such 
that e = U u V u e. Making use of e = p(l + W) and the conservation law V^(pU^) = 0(3), 
Eq. (13.1 ip can equivalently be written as 

P + f3 = -^W + 0{3). (3.12) 

Using the explicit expressions (13.9b|) . (I3.9cp . and ( 13 . 10[) for P, (3, and W, one can easily 
check that (I3.12p holds in the particular case of our dark matter model. 



IV. CURVATURE PERTURBATION UP TO SECOND ORDER 
A. Langlois-Vernizzi formalism 



We now apply the formalism introduced in Refs. [161 . Il7j for the evolution of the curvature 
perturbation. In this formalism, the usual gauge-invariant curvature perturbation of linear 
perturbation theory is generalized to a covector (or one-form) for which an exact evolution 
equation is found, valid at any perturbative order and at all scales. First, we define the local 
number of e-folds a by = 3d, and a derivative operator projected orthogonally to a given 
four- velocity field, namely _D M =^L^ U, W V . The definition of the curvature one-form is then 

C M = iV*- j£> M e. (4.1) 

Notice that one could replace the projected derivatives by covariant derivatives V M , or even 
ordinary gradients in the above formula. Using the (approximate) conservation law (13.111) . 
one obtains the following evolution equation: 

6 V P + \ 

^uC, = —(d,(P + 0)- D,e J . (4.2) 
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We refer the reader to Ref. [17| for the original derivation of this relation, but we provide 
in Appendix |A] an alternative proof based on the Cartan identity. The RHS of (14. 2 j) would 
vanish in the case of a perfect, pressureless fluid, for which f3 = P = 0, so that is 
conserved along the worldlines of CDM particles. In our case, at second order in cosmological 
perturbations, the difference between cold dark matter and dipolar dark matter is entirely 
contained in the (standard plus dissipative) pressure term (13.121) . Now, notice that the 
projected derivative picks out the spatial derivative at this order, so for instance D^e = 0(1). 
Making use of the fact that W = 0(2) already, we obtain 

&utn = \d,W + 0(3) , (4.3) 

which can then be recast, using again = 0(1), into a more familiar conservation law for 
the variable = £ M — ^D^W: 

&vl* = 0(3) . (4.4) 
One can also introduce a modified local number of e-folds, 5 = a — \W , such that the 
one- form has a definition similar to that of namely 

~ a 

= D^a — — D^e . (4.5) 

The conservation law (14.41) for the one-form defined with respect to a can be understood 
as follows: introducing ( 13.121) into ( 13. lip , one gets e + (9 — W) e = 0(3), which is the 
continuity equation of a fiducial pressureless perfect fluid of expansion = — W. This is 
precisely the expansion associated with the local number of e-folds a — the definition of O 
as the expansion of the velocity does not intervene in the derivation of the relation ( 13.121) . 



B. Curvature perturbation to second order 

We now apply the results of the calculations given in [161 for second-order cosmological 



perturbations. We use a FLRW background with line element ds 2 = a 2 (t])(— &rj 2 +^ij&x l dx : >) 1 
where a(r]) is the scale factor and 7^ the spatial metric, and we denote by ' the derivative 
with respect to the conformal time rj. Using the conservation law (14 .3p at first order, namely 

CI = 0(2) , (4.6) 

we obtain, at second order, the following expression for introducing the velocity pertur- 
bation v % as SU l = v l j a: 

Co = ~v l di( + 0(3) , (4.7a) 

O = cK + 0(3). (4.7b) 

Here ( is a scalar that generalizes the curvature perturbation scalar, of which the first and 



second order expressions in terms of the perturbations of a and e can be found in Ref. [16 
and are recalled in Appendix [B] below. Specializing the expression (jB6]) for the Lie derivative 
to second order yields, in combination with ( I4.3p . 

C + d l (v^ j ) = ^d t W' + 0(3). (4.8) 
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Then, using Eq. f!4.7bj) we get, up to an irrelevant homogeneous term, 



C + v%( = ^ + 0(3) . (4.9) 

On super-Hubble scales, the term involving a spatial gradient can be neglected, and we find 
[up to corrections 0(3)] 

W 

This shows that, on large scales, the only second-order effect of our modified dark matter 
model, compared to a standard CDM model, is to correct £ by the quantity W/3. Denoting 
by Ccdm the conserved curvature perturbation generated within the standard CDM scenario, 
we thus conclude 

W 

C = CcDM + y (4.11) 

The additional effect, where W is given by Eq. (13 . 1 0[) with (13. ip . is the sum of two terms 
which scale like the norms squared of the space-like dipole moment and its time derivative; 
schematically W ~ a£± + i±. 



C. Time evolution of the vector perturbations 

Since the internal energy W is quadratic in £^ and k m = working at second order 

in perturbation theory allows us to use the first-order results for the evolution equations, as 
previously derived in Ref. [14]]. Following their notations, we define 

e ± ={0,y), (4.12) 

with A* = 0(1) the perturbative spatial dipole moment. Easy calculations yield the following 
expressions: 

^ = I (0,\ i ' + m l ) +0(2), (4.13a) 

- (0^') +0(2), (4.13b) 
a 



k" 



^ = \(0,X i ")+O(2), (4.13c) 

where % = a' /a is the Hubble parameter in conformal time rj, and we recall that ' = d/drj. 
Using an overbar to denote background quantities, we get from Eq. ( 13 . lOj) : 

W = 2<iraa 2 \ i \ i - -A i; Aj + 0(3) , (4.14a) 

W = — (A l, A' - 2naa 2 \ l \i) + 0(3) , (4.14b) 
a ' 



where we denote \ = 7jjA J and A^ = 7ijA J/ . In order to derive Eq. (I4.14bj) . we used the 



first-order accurate evolution equation for the dipole moment X\ which takes the remarkably 



simple form given in Eq. (3.29) of Ref. [14 



A 4 " + HX 1 - 4iraa 2 \ l = 0(2) . (4.15) 
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This equation can be solved, at first order, using the FLRW equations for the background 
dynamics. If we consider the background to be successively dominated by radiation and by 
(dipolar) cold dark matter, the scalings for the background energy density are p r oc a -4 for 
the radiation domination era and a oc a~ 3 for the matter domination era, which allows us 
to find the approximate solutions 

A* (77, x) = ^(x) J ( v /6O77) + B*(x) K (^h]) (radiation domination) , (4.16a) 

A l (?7, x) = C"(x) rj 2 H j-^- (matter domination) . (4.16b) 

Here Iq and K are modified Bessel functions of the first and second kind, respectively, and 
A\ B\ C\ and D % are (space-dependent) integration constants. We introduced the inverse 
of a characteristic timescale, £1 = a eq (87rp e q/3) 1 / 2 , the subscript "eq" indicating the time 
when the background densities of the radiation and matter fluids are equal, and p eq being 
the common value of the radiation density p r and the matter density a at equilibrium. 

In fact, it is possible to give a fully analytic solution to the case of the mixture of the two 
fluids. For this purpose, we rewrite the dynamics in terms of the variable y = a/a cq . The 
FLRW equations for the background are then given by 

H 2 = fi2 y + 1 ; (417a) 

ft = 2~> 4.17b 

2 y 2 

where we used the equation of state P r = p r /3 for the radiation. Finally, we get an evolution 
equation valid in the two regimes, radiation as well as matter eras [hereafter we drop the 
neglected remainder 0(2)]: 

2y(y + 1)_ + (3y + 2) — - 3A l = . (4.18) 

This homogeneous, second-order differential equation is called Meszaros equation. It appears 
in standard cosmology in the context of the growth of the matter density contrast on sub- 
Hubble scales; see e.g. Ref. [2l|. Its most general solution is known analytically, and reads 

A%, x) = A*(x) X + (y) + ^(x) \~(y) , (4.19) 

where A 1 and B l are again (space-dependent) integration constants, while the growing (+) 
and decaying (— ) modes read 

A + = y + I , (4.20a) 



A" 



VTFi-(. + !)>.($g±i). (4.20b) 



One can check that Eqs. f!4.20p reproduce the asymptotic behaviors of the solutions (j4. 16aj) 
and ( 14.16bj) . In particular, \~(y) behaves like ln(y) [or ln(r^)] when y — >■ 0, and like y~ 3 ^ 2 



(or r/~ 3 ) when y — > +00. We thus have an analytic expression for X(y), which we can inject 
into the formula fl4.14a[) for W. This immediately gives, via Eq. (14. lip , the time evolution 
of the curvature perturbation ( on super- Hubble scales. 
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Discarding the decaying mode A , we can write the evolution of the dipole as 

A*(y,x) = (l + |y) Al(x), (4.21) 

where A* (x) is a time- independent vector field corresponding to the value of the dipole long 
before the equivalence matter-radiation (j/ < 1). Substituting into the expression (|4.14a|) 
for W, one obtains the time evolution of W as a function of the "primordial" dipole A* : 

W = J^(a)Xl, (4.22) 

where A 2 = 7^ A* X{, and the evolution factor is given in terms of the variable y = a/a eq by 

This quantity diverges when going backwards in time, and we will use it only sufficiently 
late so that W remains small enough to be treated as a second order perturbation of Ccdm- 
After equivalence, J 7 will increase linearly with the scale factor. At the radiation-matter 
equality, the typical value of the correction W is given by ~ (% eq A*) 2 , which is small if the 
primordial dipole moment is small compared to the Hubble radius at equivalence. 



V. STATISTICAL DESCRIPTION OF THE VECTOR PERTURBATIONS 
A. Formalism and notations 

In this section, we examine how statistical fluctuations of the DM dipole could affect the 
nonlinear cosmological perturbations and therefore leave a potentially observable signature. 
To do so, we will use the extension (22L 23 1 of the SN formalism j24-23|, where the curvature 



perturbation ( is allowed to depend on a vector field Aj. In the present case, we decompose 
the total curvature perturbation ( into the standard CDM contribution and a contribution 
due to a vector field A^ directly related to the dark matter dipole, so that one can write, 
up to second order, 

C = Ccdm + N* 5A t + ^ 6 A SA d , (5.1) 

where we use implicit summation over repeated indices. 

Let us now characterize the statistical properties of the vector perturbations, following 
closely the general approach developed in Refs. [22I, 23[. Since the space-like vector field 



introduced in the present model is a priori purely phenomenological, we will simply assume 
that it is a collection of uncorrelated Gaussian random fields, without trying to determine its 
spectrum from a more fundamental description. 4 It is sometimes convenient to decompose 



4 Such a description would consist, for instance, in embedding the dipole in an inflationary model and 
computing its spectrum from the amplification of its quantum fluctuations during inflation, as was done 
for fundamental vector fields in Ref. [28[. 
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the three components of the vector field, in Fourier space, into the familiar longitudinal (£), 
left (L), and right (P) polarizations, according to the expression 



6M V ,k)= ( d 3 xe- ikx ^(r / ,x) = ^M a (r ? ,k)e°(k), (5.2) 

where the polarization vectors read ef = (1, i, 0) / \/2, ef = (1, — i, 0)/v2, and e\ = (0,0,1), 
with the z direction aligned with the direction k of the wavevector k. Following the conven- 
tions of [22], we let the change k — > — k reverse the z and x directions, but not the y direction, 
such that ef*(k) = — ef (— k). Together with the reality condition 5A*(k) = 5Ai(— k), this 
implies 5A* a (k) = -5A a (-k). 

The random fields 5A a (k) are assumed to be Gaussian and statistically isotropic, with no 
correlation between different polarization modes or with other fundamental fields (such as 
the inflaton, the curvaton, etc.). The power spectra P a (k) for each polarization are defined 
by the relations 

(SA a (k) 5A p {k!)) = -(2vr) 3 5(k + k') P a {k) 8 aP , (5.3) 

where 5 is the usual three-dimensional Dirac delta function, 5 a p is the Kronecker symbol 
and k — |k|. With these conventions, the power spectra P a (k) are positive-definite. Using 
the decomposition (15.21) . the two-point correlation function of the dipole can be written, in 
Fourier space, as |22| 

(SMk) SAjik')) = (27r) 3 (5(k + k / )P i ,(k) 

= (2vr) 3 5(k + k')[7- vcn (k)P + (A;) + i P° dd (k)P_(A;) + T^(k)P e (k)] , (5.4) 

with the even-parity, odd-parity, and longitudinal tensors defined by T^ ven = Sy — k{kj, 

T° dd = Eijkkk, and T^ ong = kikj, respectively. In the previous expressions, we introduced 
the usual totally antisymmetric Levi-Civita symbol Syjj., such that £123 = +1, as well as the 
power spectra 

P ± = l -{P R ±P L ). (5.5) 

Note that the Lagrangian (I2.2p of DDM does not break the parity symmetry. If this property 
is also shared by the Lagrangian describing the dipole at a more fundamental level, then the 
left (L) and right (P) power spectra Pl and Pr should coincide, implying P = 0. In the 
following, we shall keep a general P_. 

In the rest of this section, we will apply the above formalism to the DM dipole and 
compute the power spectrum and bispectrum of the curvature perturbation. Since we have 
considered the dipole as a perturbation with vanishing background value, we have a priori 
N 1, = 0, which implies that the tree-level contributions to the spectrum and bispectrum 
simply vanish. One must then consider the one-loop contributions, as we do in the next 
subsection IV Bl However, in the subsequent subsection IV C\ we will instead assume that 
the dipole can be decomposed into a small "homogeneous" component, corresponding to its 
average value in our accessible region, and a spatially-dependent component, in which case 
the term linear in 5Ai in Eq. 05. ip is nonzero. We will check that the tree-level expressions 
for the spectrum and bispectrum are compatible with the one-loop calculations of Sec. IV B\ 
in the sense of cosmic variance. 
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B. One-loop contributions to the spectrum and bispectrum 

Comparing Eq. ( 15. ip with the expressions (I4.1ip and (I4.22p . one can directly identify the 
vector degree of freedom 5Ai with the primordial DM dipole A* , which we have treated as a 
perturbation, so that 

5A = X„, N l = 0, N^ = ^(a)5 l \ (5.6) 

where we assume 7^ = 5ij (flat spatial sections) for simplicity. At the tree-level, the total 
power spectrum is given by 

P f (k) = P CcDM (k) + AWP^k) , where (C(k)C(k')) = (2vr) 3 5(k + k')P c (k) . (5.7) 

Because N l — the contribution of the vector field vanishes. The first non-vanishing effect 
of the dipole appears in the one-loop contribution 

Pr l0 ° P ( k ) = \ N ' l3Nkl j (^3 ^(P)^(k - P) . (5.8) 

As an illustration, we will assume in the following that the power spectrum for the dipole is 
scale invariant, i.e. Py(k) oc k~ 3 . Then, the above integral is divergent at the two singular 
points p = and p = k. These logarithmic divergences can be regularized by introducing 
a cut-off L^ 1 , where L is a length scale which, in practice, can be chosen as the present 
Hubble radius: L = Hq 1 . For instance, the singularity at p = gives 

r k d 3 p 4tt 

/ -f Vijip) = — (2P + + Vt) ln(fcL) Sij , (5.9) 

where we denote Pij(k) = Py(k) k 3 /(2iT 2 ) and similarly for P +j ^(k). This gives the following 
estimate for the one-loop contribution to the power spectrum: 

^-^(k) ~ ^rF 2 (2V + + V f fk~ 3 ln(kL), (5.10) 

where we have substituted N l i = |j r 5 i - 7 . This must remain smaller than P^ CDM in order to 
satisfy observational constraints. 

We now turn to non-Gaussianities. The bispectrum of the curvature perturbation, 
P^(k, k', k"), corresponds to the three-point correlation function in Fourier space, 

(C(k)C(k')C(k")) = (2nfS(k + k' + k")B c (k, k', k") . (5.11) 

Similarly to the power spectrum, the tree- level bispectrum vanishes when N l = 0, whereas 
the one-loop contribution is given by 

^- loop (k, k', k") = N^N kl N mn J ^ P fe (p)P jm (k - P )P n (k' + p) . (5.12) 

With a scale-invariant spectrum Vij, the above integral contains logarithmic divergences for 
the values p = 0, p = k and p = — k'. Introducing as before the cut-off L~ l , one finds for 
the bispectrum the approximate expression (valid at leading order when L — > +00) 

flj-^k, k', k") ~ N^N hl N mn (2P+ + V e ) HkL) 5 ik P jm {k)P ln {k') 

= |-^ 3 (2P+ + Vt) HkL) P tJ (k)P l3 (k>) . (5.13) 
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C. Tree-level spectrum and bispectrum 



An alternative derivation of the spectrum and bispectrum consists in assuming that the 
dipole A*, which we have considered as purely perturbative so far, can be decomposed into 
a nonzero homogeneous component and a spatially-dependent component. Physically, the 
homogeneous component can be interpreted as the background vector in our accessible box 
of size Hq 1 and corresponds to the average field over our observable region in the particular 
realization associated with our Universe. It has to be small in order to be consistent with 
the analysis of Sees. [TTT] and HV] and to be compatible with current bounds on anisotropy in 
the CMB. Substituting the decomposition 

a;(x) = a: + 5a;(x) (5.i4) 

into the formula (I4.22p . the curvature perturbation ( 14.1 ip can be written in the form 

C = Ccdm + + S\» + - . (5.15) 

The second term is an irrelevant homogeneous contribution that can be discarded, and the 
above expression is thus of the form (15. ip with a nonzero linear term: 

6Ai = 5\*i , N* = , N** = \T0 . (5.16) 

Assuming that Ccdm and 5Ai are uncorrelated, the tree-level power spectrum for the curva- 
ture perturbation is given by 5 

P c tree (k) = P CcDU (k) + iV'AP [T°™(k)P + (k) + T^(k)P £ (k)] 

= Pf° (k)[l + g (N%) 2 ] , (5.17) 

where the angle-averaged, isotropic spectrum P™°(k) and the scale-dependent, anisotropy 
fraction g(k) read (posing N 2 = N l N l ) 

PT = P <c D m+N 2 P + , (5.18a) 
9 = N 2 ^-^. (5.18b) 

Present observations of the CMB temperature fluctuations require g to be small enough, 
typically g < 0.1. 6 In our case, consistency with our assumption that the dipole modifies 
cosmological perturbations only at second order implies that g is small, such that P^ ~ P™°, 
and that the isotropic spectrum P™° is dominated by the standard Ccdm contribution, namely 
piso r^, P^ cdm . Note that in the particular case where Pi = P + , the power spectrum is 
automatically isotropic. 



5 The term involving P_ disappears because N l N J is symmetric. 

6 Although an anisotropy has been detected in the WMAP data, the direction of this asymmetry coincides 
with the direction orthogonal to the ecliptic, suggesting that this is a systematic effect; see Refs. 



29, 



14 



Let us now turn to the bispectrum of the curvature perturbation. Since the coefficient N l 
is now nonzero, the tree-level contribution no longer vanishes, and is given by the expression 

B? ee (k, k', k") = N' l WN kl [P lfc (k)P j/ (k / ) + 2 perm.] 

= KX [P ik (\c)P jk (k') + 2 perm.] . (5.19) 

One can show that this result agrees with the 1-loop contribution of the previous subsection 
IV Bl by recalling that N l = f-^A* corresponds to the background vector in our accessible 
box of size Hq 1 . It can be interpreted as the average dipole field over our observable region. 
One can estimate the typical value of this background field — or rather of the tensor A* A{ 
which appears in the bispectrum — by taking the variance matrix 

<(A:(x)Ai(x)» L = ^ P^(p) = l - {2V + + V t ) ln(fc max L) , (5.20) 

where fc"^. corresponds to the smallest accessible length scale. Substituting A* A^ in the 
bispectrum (15 . 19[) by the above expression, one recovers the one-loop bispectrum ( 15 . 13[) at 
the leading order when L — > +oo. 7 This shows that the computation based on a non-zero 
background value for the dipole moment is also compatible, in the sense of cosmic variance, 
with a zero averaged value for the background dipole moment. 

The shape of the bispectrum (I5.19P naturally involves the preferred cosmological direction 
A*. In general it is highly anisotropic and quite complicated. In the particular case where 
Pe = P + and P_ = 0, the bispectrum becomes isotropic with local shape: 

P c (k, k', k") = - X 2 [P+{k )P + (k') + P + (k')P + (k") + P + (k")P + (k)\ . (5.21) 

It is then entirely characterized by the isotropic spectrum P + (k), grows like J 73 (a), and is 
proportional to the norm A^ of the primordial homogeneous dipole. 



D. Discussion and conclusions 

In order to get an estimate of the amplitude for non-Gaussianities, it is useful to introduce 
the usual parameter 8 

/NLl ' ' ' ~ 6 Pc(k)P c (k') + P c (k')P f (k") + Pc(k")P c (k) ' {b - ZZ) 

In general /nl has a complicated (anisotropic) structure in Fourier space. To determine its 
order of magnitude, let us consider the special case Pi = P + , where the bispectrum shape is 
of the local form (I5.2ip . and let us assume that P+{k) is quasi-scale- invariant, like P( CDM (k). 



7 Similarly one can also recover by the same calculation the spectrum (|5.10p . 

8 The factor 5/6 in Eq. (|5.22[) comes from the original, historical definition of /nl as a parameter measuring 
non-Gaussianity of the local type in the Bardeen gravitational potential $ = |C, according to the formula 
$(x) = $g(x) + /nl ($g(x) - ($g( x )))' where $g(x) is a Gaussian random field [HlJ - 
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Then P((k) is also quasi-scale-invariant, and we have for the non-Gaussianity induced by 
the dark matter dipole 

f^lTni P f r %, (5,3) 

where W = J^Xl is the background value of the correction (I4.22p to the CDM curvature 
perturbation, and we introduced the fraction S of the curvature spectrum due to the dipole, 
namely [remember Eq. (15. 18ei[) 

E^N^ = ^Xl^±. (5.24) 

The dimensionless factor S is assumed to be very small in the present work, and thus tends 
to suppress the amplitude of non-Gaussianities. However, the denominator W in Eq. (15. 23ft 
is also assumed to be very small (W <C 1CT 5 ) in order to recover the standard results of linear 
cosmological perturbations. Therefore, the amplitude of the non-Gaussianities produced by 
the dipole moment could be quite significant. Furthermore, we find that this amplitude 
increases with time after the radiation-matter equality, roughly like J 73 ~ a 3 , because the 
dark matter dipole itself grows after that epoch. This behavior is in sharp contrast with 
usual models of primordial non-Gaussianities, where ( remains frozen on super-horizon scales 
during the standard cosmological era. 

That untypical prediction could be used to test the possible dipolar nature of dark matter 
as follows. Non-Gaussianity in the curvature perturbation can be measured independently 
through the bispectrum of the temperature fluctuations in the CMB, as probed by the Planck 
satellite, and the bispectrum of the mass density distribution, at low redshifts (z < 2) and 



large scales, as probed by future galaxy surveys (e.g. Euclid, SDSS-III) |32|, [33|. In order 
to make specific predictions for the relative amplitudes of the non-Gaussianity that could 
be measured by these various probes, one should extend the analysis of the present work by 
taking into account the presence of the dipole in the calculation of the transfer functions. 

To summarize, we investigated the model of DDM at the level of second-order cosmolog- 
ical perturbations. We showed that, on large scales, the only effect of the dipole moment is 
to correct the curvature perturbation generated within the standard scenario, say Ccdm, by 
a term proportional to the internal energy W of the dipolar fluid, which is quadratic in the 
dipole moment [see Eq. (14. lip ]. We then derived, by two different methods in Sees. IV Bl and 
IV CI the amplitude of the non-gaussianities generated by this specific contribution. Unfor- 
tunately, without a fundamental understanding of the origin of the dipole field introduced in 
this phenomenological model of dark matter, the amplitude of the dipole moment at early 
cosmological times cannot be determined, leaving the global scaling of the additional signal of 
non-Gaussianity (1 5.231) — (15. 241) unconstrained. However, the property that the (second-order) 
contribution of the dipole to the curvature perturbation is time-evolving on super-Hubble 
scales would lead to specific relations between the levels of non-Gaussianities observed by 
various probes. 
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Appendix A: Alternative derivation of the relation (|4,2p 

We provide here an alternative and more compact derivation of the Langlois-Vernizzi re- 
lation (14. 2p . based on the use of differential forms and the Cartan identity. In the following, 
boldface letters indicate differential forms, arrows indicate four-vectors, and the dot stands 
for the contraction of a four-vector with the first index of a differential form, while d repre- 
sents the exterior derivative operator. Replacing the projected derivatives D M by ordinary 
derivatives in the definition (14. ip of and using Eq. (13. lip to rewrite the ratio a/s, we 
may express Eq. (14. ip as 

£ = da + - — -. (Al) 

S 3(e + P + (3) K } 

We will make use of Cartan's identity for the Lie derivative of a differential form, which 
reads for a four- vector V and a differential form u> of any rank, 9 

^u = V -du> + d(V -u). (A2) 

We are going to apply this identity to the four-vector U and the one-form £. We first notice 
that, by construction (see Eq. ( 14. ip again), £ is orthogonal to the four- velocity: U • C, = 0. 
We also note that, in the particular case of a pressureless perfect fluid with P = (3 = 0, £ 
becomes an exact one-form, £ = d(a + | hie), and we readily obtain that the RHS vanishes. 
Coming back to the general case, we have 

d < = d ( da+ 3( £ + P +fl ) dg )— g d(P + ^ )Ad£ - (A3) 

where we used Eq. (13. lip again. Now, contracting with U, we obtain the desired relation: 

V P + \ 

Sf ff C = gj (d(P + 13)- — p de) , (A4) 

which is identical to Eq. (14. 2p . 



Appendix B: Expressions for second-order perturbations 

For completeness, we provide here some explicit expressions for the curvature perturba- 
tion in terms of first- and second-order quantities in cosmological perturbations. They can 
be found in Ref. [l6j], together with other definitions and explicit expressions for the metric 
perturbations. We adopt the convention that any quantity Q is to be decomposed in the 
following way: 

Q( V , X ) = Q( V ) + Q {1) ( V , X ) + Q i2) ( V ,^) + 0(3), (Bl) 



9 A consequence of this identity is the formula d(Jz?^ uj) = Jz?^ (du>). 
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where Q stands for the background quantity and Q(i), Q(2) for the first- and second-order 
perturbations, respectively. Notice the change in convention with respect to Ref. 16] . where 
the authors keep a factor 1 /2 in front of the second-order perturbation in Eq. (1B1I) . 

For a FLRW background, we have ds 2 = a 2 {rj) {—drf + 7^ dx l dx 3 ') . We recall that a 
is the integrated expansion, defined by 6 = 3d (see Ref. [l6[ for an explicit expression in 
terms of metric and matter perturbations), and that we define Ufo = v l /a\ The curvature 
one-form vanishes identically in the background, and we find at first order 



Co 



(i) 



c 



(1) 



0, 



(B2a) 
(B2b) 



with the following definition for the scalar curvature perturbation (: 



1) 



a 

"(l) - -7^(1) 



(B3) 



where the symbol ' stands for a derivative with respect to the conformal time 77. We see that 
at first order (0 vanishes and that Q can be written as a pure spatial gradient. This is not 
the case in general, and at second order we get 



.(2) 
.0 

.(2) 



(B4a) 
(B4b) 



with 



a' 



1 , a! , 1 

C(2) = «(2) - — £(2) - pa(i)£(i) + |^ £ (i) £ (i) + 7^7 



— e 



(i) 



(B5) 



On large scales, these expressions agree with other expressions used in the litterature for 
the conserved curvature perturbation; see Ref. 16] for a discussion. We also reproduce here 
the expression of the Lie derivative of Q along at second order, introducing A as the 
(first-order) perturbation in the goo component of the metric as goo — — a 2 (l + 2 A): 



(2) 



Si 



(2)' 



(B6) 



Since C| is a pure gradient, the last two terms in the RHS of Eq. flB6j) can be combined to 
give di(y J Q )> as used in the evolution equation f )4.8p . 
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